Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 4 March 2009 (MN KTeX style file v2.2) 



Segue 2: A Prototype of the Population of Satellites of Satellites 



V. Belokurov 1 , M.G. Walker 1 , N.W. Evans 1 , G. Gilmore 1 , M.J. Irwin 1 , M. Mateo 2 , 
L. Mayer 3 , E. Olszewski 4 , J. Bechtold 4 , T. Pickering 5 . 



On 
O 
O 
(N 



< 
O 

6 



> 

00 
OC 

q 
o 

On 
O 



- 1—1 

x 



1 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge, CB3 0HA, UK 

2 Department of Astronomy, University of Michigan, Ann Arbor, MI 48109, USA 

3 Institute for Theoretical Physics, University of Zurich, CH-8057 Zurich, Switzerland 
^Steward Observatory, University of Arizona, Tucson, AZ 85721, USA 
5, 



MMT Observatory, University of Arizona, Tucson, USA 



4 March 2009 



ABSTRACT 

We announce the discovery of a new Milky Way satellite Segue 2 found in the data of Sloan 
Extension for Galactic Understanding and Exploration (SEGUE). We followed this up with 
deeper imaging and spectroscopy on the Multiple Mirror Telescope (MMT). From this, we 
derive a luminosity of M v = -2.5, a half-light radius of 34 pc and a systemic velocity of 
~ -40 kms -1 . Our data also provides evidence for a stream around Segue 2 at a similar 
heliocentric velocity, and the SEGUE data show that it is also present in neighbouring fields. 
We resolve the velocity dispersion of Segue 2 as 3.4 kms -1 and the possible stream as ~ 7 
kms -1 . This object shows points of comparison with other recent discoveries, Segue 1, Boo 
II and Coma. We speculate that all four objects may be representatives of a population of 
satellites of satellites - survivors of accretion events that destroyed their larger but less dense 
parents. They are likely to have formed at redshifts z > 10 and are good candidates for fossils 
of the reionization epoch. 
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I INTRODUCTION 

The idea that the outer parts of Galactic haloes are built up from 
the merging and accretion of satellites is now well-established. The 
building blocks that contributed most to the Galactic halo have been 
broken down into streams of debris. Reconstructing the history 
is difficult as the progenitors have been disassembled and phase 
mixed. In Cold Dar k Matter cosmogonie s, smaller haloes form ear- 
lier and are denser l INavarro et alJI 19970 . So, the entourage of the 
accreted progenitors, smaller satellites of the bigge r satellites, may 
have survived against tidal destruction (see e.g., iDiemand et al.l 
l2008h . 

Amongst the recent discoveries of Milky Way satellites, 
there are objects whose properties are unlike conventional glob- 
ul ar clusters or dwarf ga l axies, such as Co ma, Segue 1 and Boo 

II teelokurov et"ai]|2007l ; Iwalsh et alj|2007h . They have half-light 
radii of ~ 30-70 pc and luminosities below M r = -3. In fact, Coma, 
Segue 1 and Boo II lie projected on the Sagittarius stream, and 
have velocities consistent with Stream membership. Irrespective 
of whether they are dwarf galaxies or globular clusters, it seems 
reasonable to conclude that they once belonged to the Sagittarius 
galaxy. Could they be the first examples of a population of satel- 
lites of satellites? 

In this Letter, we report the discovery of a further object anal- 
ogous to Coma, Segue 1 and Boo II. Using the recently available 



SEGUE imaging, we have extended our ongoing survey of stel- 
lar overdensities in the outer Milky Way halo. We followed up the 
new object - called Segue 2 - with deep imaging and high resolu- 
tion spectroscopy and present its properties and possible genealogy 
here. 



2 DATA AND DISCOVERY 

The original Sloan Digital Sky Survey (SDSS) imaged most 
of the North Galacti c Cap plus three st ripes of data in the 
Sout h Galactic Cap iAbazaiian et al.l l2009h . The SEGUE sur- 
vey jYannvetalJl2009l) " is primarily spectroscopic, but comple- 
ments the SDSS imaging data with 15 2.5° wide stripes along 
constant Galactic longitude, spaced by approximately 20° around 
the sky. The stripes probe the Galaxy at a variety of longi- 
tudes, sampling the changing relative densities of thin disk, thick 
disk, and halo. The SEGUE imaging footprint is illustrated at 
http : //www . sdss . org/dr7/seguephoto_big . gif 

Figu re Q] shows the fami liar picture of the "Field of 
Streams" teelokurov et alj2 006a). together with the SEGUE imag- 
ing stripes with right ascension a < 90° and declination 6 < 50°. 
We have also cut on Galactic latitude b < -20° to avoid showing 
regions dominated by the Galactic disk. As is usual, the magni- 
tudes of stars have been corrected for extinction using the maps 
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Figure 1. Location of Segue 2 marked by asterisk with respect to "the Field of Streams", in SEGUE imaging. This is a stellar densi t y plot o f all stars with 
20.0 < i < 22.5 and -1 < g — i < 0.6. The magnitude range is divided into three equal-sized bins analogous to Belokur ov et alj j2006al) . Note that the 
Sagittarius Stream trailing arm is clearly visible crossing the equator at right ascensions a ~ 40° . 
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Figure 2. Leftmost two panels: SEGUE image of the 5' X 5' field centered on Segue 2 and density of all objects classified as stars in a 30' X 30' field. Note 
the bright saturated star and associated artefacts. The annuli are used in the construction of the CMDs. Rightmost three panels: CMD of the inner and outer 
regions, together with their Hess difference. There is a clear main sequence turn-off, together with sparse red giant branch (RGB) and blue horizontal branch 
(BHB). 



of lSchlegeletai] ( ll998l) . Ru nning the algorithm for i dentification 
of overdensities described in lBelokurov et all {2006b), we isolate 
high significance peaks. The location of the highest is marked by a 
white asterisk in Figure [T] 

Figure[2] shows the SEGUE view of the overdensity. The left- 
most panel is a cut-out of the sky with the overdensity at the centre. 
The only discernible structure is in fact an unrelated saturated fore- 
ground star and accompanying artefacts. In the next panel, the den- 
sity of resolved stars is shown, and now the re is an evident overden- 
sity with significance ( Kop osov et al.l20 08) S = 4.7 at the centre of 
the image. Annuli are marked which are used to select stars within 
8' and outside 12' for the two color magnitude diagrams (CMDs). 
The rightmost panel gives the Hess difference, and there is a clear 
main-sequence turn-off, together with definite hints of a red giant 
branch (RGB) and a blue horizontal branch (BHB). 



3 FOLLOW-UP 
3.1 Imaging 

Follow-up imaging of Segue 2 was carried out on Oct 7 2007 
using the Megacam imager (Mcleod et al 2006) on the Multi- 
ple Mirror Telescope (MMT). Megacam comprises 36 2048x4608 
E2V CCDS. With 2x2 binning, pixels are 0.16" and each im- 
age is 24' x 24' in size. Six 300s exposure images in r' and 7 
in g' were collected with dithers of 100-200 pixels in each co- 
ordinate between frames. Frames were processed using SAO's 
MEGARE D package and comb ined using SWARP package in 
TERAPIX ( Radovic heTal.120011) . A final set of object catalogs was 
generated from the stacked images and objects were morphologi- 
cally classified as stellar or non-stellar (or noise-like). The detected 
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Coordinates (J2000) 
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v Q , <x 

Position Angle 

Ellipticity 

rj, (Exponential) 
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182° ± 17" 
0.15 ±0.1 
3f4±0;2 
17 m 7±0 n, l 
-2?5 ± 0?2 
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Table 1. Properties of the Segue 2 Satellite 



objects in each passband were then merged by positional coinci- 
dence (within 1") to form a combined g, r catalogue and photomet- 
rically calibrated on the SDSS system using stars in common. 

Figure [3] compares the original SEGUE data with the Mega- 
cam follow-up. From the CMDs, it is immediately clear that the 
Megacam data probe at least 2 magnitudes deeper and reveal t he 
main sequence of Segue 2 . The algorithm o flMartin et al] ( 120081) is 
applied to stars selected by the masks shown in the Figure. Model 
isodensity contours, based on both the SEGUE and Megacam data, 
are shown in the lower panels. We see that the shallower SEGUE 
data yields a slightly more extended and elliptical light profile. 
The extracted structural parameters using the deeper Megacam data 
are listed in Table [T] There are possibly 4 BHB stars associated 
with Segue 2 , which can be used to obtain a distance modulus 
m — M = 17.7 or 35 kpc. Using this, the half-light radius is r h = 34 
pc which is small compared to typical ul trafaint dwarf galaxies 
( iBelokurov et alj|2 007; Gilm ore et alj2007h . 
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Figure 4. Left: CMD of all stars (black dots) within 2r^, with objects followed-up spectroscopically circled. Filled circles with valid velocities and ZMg color- 
coded such that gray shows the foreground population, red shows red clump giants, blue show BHBs and orange subgaints in Segue 2 . The M92 ridgeline 
(including the BHB) is overplotted together with the M13 ridgeline for comparison. Middle: The distribution of ZMg values, showing clean separation between 
dwarfs in the thick disk (high XMg ) and giants in Segue 2 (low SMg ) at bright magnitudes. The vertical solid line marks the boundary adopted here. Note 
that at fainter magnitudes, there is a population of stars with intermediate ZMg and velocities similar to the systemic velocity of Segue 2 . Right: Heliocentric 
radial velocities of all stars with good spectra. The characteristic velocities associated with the halo and thick disk at this location are enclosed by vertical 
dashed lines, together with solid lines showing the systemic velocity of Segue 2 plus or minus 10 kms~'. 
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Figure 3. Top left and right; CMDs of Segue 2 in SEGUE and in the follow- 
up Megacam data, which goes at least 2 mag deeper. The follow-up data 
shows a well-defined main sequence, with grey lines marking the ridgeline 
of M92 and the mask used to select members. Bottom left and right: Grey- 
scale density of stars selected by in the 24' X 24' field of view centered on 
Segue 2 . The ellipses show the model isodensity contours corresponding to 
rt, and 2r^. The colored dots show the locations of likely members followed 
up spectroscopically. 



3.2 Spectroscopy 

On 22-23 and 26 October 2008, we obtained high resolution spectra 
of 358 targets around Segue 2 using three independent fiber config- 
urations with the Hectochelle spectrograph at the MMT. Over the 1° 
field, we targeted RGB candidates as faint as i ~ 21.5, as well as the 
handful of BHB candidates at ( ~ 18.5. We ext racted and ca l ibrate d 
spectra following the procedure described by iMateo et al.l (2008). 
The Hectochelle spectra span 5 150-5300 A, where RGBs and fore- 
ground dwarfs prominently exhibit the Mg-I/Mg-b triplet (MgT) 
absorption feature. For the RGB candidates, we measure velocity 
by cross-correlating each spectrum against a high signal-to-noise 
template spectrum, built from co-added spectra of late-type radial 



velocity standards observed with Hectochelle. However, this tem- 
plate poorly resembles the spectrum of a BHB, in which high tem- 
perature suppresses Mg absorption and the only prominent absorp- 
tion feature is the Fel/Fell blend at 5169 A. Thus, for each BHB 
candidate, we measure the centroid of this feature and calculate ve- 
locity directly from the redshift. For all spectra, we also measure a 
composite magnesiu m index, £Mg, effec tively a pseudo-equivalent 
width for the MgT ^Walker et al]|2007l) . This quantity correlates 
with metallicity, temperature and surface gravity, and helps to sepa- 
rate members from foreground. We determine measurement errors 
fo r both velocity and £Mg using the bootstrap method described 
bv lWalker etaIlj2009h . Velocities and EMg in the Segue 2 sample 
have mean (median) errors of 1.1 (0.6) km s _I and 0.20 (0.17) A , 
respectively. 

Figure [4] shows correlations between the photometric and 
spectroscopic properties of stars within 2r h . The first panel shows 
the locations of stars targeted for follow-up on the CMD, together 
with the ri dgelines of M92 ([ Fe/H] = -2.24) and M13 ([Fe/H] = 
-1.65) from IClem et al.l j2008h . The second and third panels show 
EMg and i'hel plotted against ( band magnitude, from which we 
will identify three classes of Segue 2 members, namely bright red 
clump giants (RCGs), fainter RGBs and BHBs, and two types of 
contaminants, Galactic foreground and possible tidal stream. 

One population clearly stands out in the EMg panel. Given 
the magnitude distributions, it is evident that the stars with higher 
SMg are dwarfs in the Galactic thick disk. At a similar magnitude 
i ~ 17.5, there is only one other population with a tight EMg distri- 
bution, namely the red clump stars visible in the first panel. These 
form a narrow distribution in velocity space, suggesting that the 
systemic velocity is ~ -40 kms~'. Sample spectra of a Galactic 
dwarf and a RCG are shown in Figure [5] Note that the dwarf has 
higher surface gravity and potentially higher metallicity than the 
RCGs, and hence has broader absorption features. The velocity sig- 
nature is corroborated by the three BHBs, which all have velocities 

40 kms~' . Finally, we use the systemic velocity of Segue 2 with 

a range of ±10 kms~' as a secondary cut, as shown by the vertical 
lines in the third panel. This gives a further three candidate mem- 
bers (shown in yellow), which all lie redwards of the M92 ridgeline 
and are better described by more metal-rich templates like M13. An 
example of their spectra is illustrated in the fourth panel of Figure[5] 
Unfortunately, the low signal-to-noise does not allow the direct ex- 
traction of the metallicity for these stars. 
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Figure 5. Examples of spectra of thick disk contaminant (top), together 
with three members of Segue 2 , a RCG (upper middle), a BHB (lower 
middle) and a faint RGB (bottom). Only the bluer part of the wavelength 
range containing the MgT is shown. 
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Figure 6. Heliocentric radial velocity versus elliptical radius (left) and SMg 
(right). All the stars with measured XMg are shown as black dots, whilst 
stars satisfying ZMg < 0.7 are circled. The extent of Segue 2 is illustrated 
by the vertical lines showing t\ and 2l\. The open circled stars do not ap- 
pear to be kinematically associated with Segue 2 . The structural parame- 
ters inferred from the photometry are consistent with the spectroscopic and 
kinematic signal. Note that the velocity histogram on the extreme right re- 
veals a substantial population of stars with -60 < vhel < - 30 kms that 
cannot be attributed to either thick disk or halo. 



These three fainter giant stars may be representatives of a dis- 
tinct population with low to intermediate IMg , clearly visible in 
the second panel. The velocities of the stars in this population are 
offset from the thick disk, but roughly centered on Segue 2 . The 
population is much more evident in Figure(6] which shows all stars 
with measured velocities and EMg within 45'. The velocity his- 
togram along the rightmost vertical of the plot shows that the fore- 
ground consists of the expected thick disk and halo populations, 
together with an unknown component, perhaps a tidal stream. The 
stars in this component extend over the entire field and they are 
kinematically colder than the thick disk but significantly hotter than 
Segue 2 . 

The left panel of Figure [7] shows the SEGUE spectroscopic 
footprint around Segue 2 . The right panel shows the distribu- 
tions of radial velocities corrected to the Galactic Standard of 
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Figure 7. Left: SEGUE spectroscopic fooprint in the vicinity of Segue 2 
. Right: Histograms of radial velocities corrected to the Galactic Standard 
of Rest in fields coded according to color. The black vertical lines shows 
the systemic velocity of Segue 2 , whilst the colored ones show the mean 
velocity of the thick disk. The smooth curves show the sums of Gaussian 
velocity distributions of thick disk «v) = 180, <x = 50 kms -1 ) and halo 
((v) = 0, <x = 100 kms -1 ). Note the prominent bumps in the red and orange 
histograms corresponding to the Sagittarius trailing arm, and similarly the 
bumps in the pale blue and green histograms. 

Rest are shown for each color-coded field. The underlying smooth 
curves come from a simple model of the Galactic thick disk and 
halo represented as two Gaussians with means and dispersions of 
«v> = 180, tr = 50 kms -1 ) and «v> = 0,cr = 100 kms -1 ) re- 
spectively. We immediately notice the presence of bumps or cold 
features in several fields. At low declination, the Sagittarius trail- 
ing arm is detected at vqsr ~ -120 kms -1 . At higher declinations, 
there are significant deviations from the smooth curves at velocities 
Vgsr ~ 40 kms -1 , which is close to the systemic velocity of Segue 
2 . 

More speculatively, we suggest that the stream-like overden- 
sity seen in the MMT data is part of the same structure as the 
features seen in the velocity histograms in the higher declina- 
tion (lower Galactic latitude) SEGUE fields. This could be con- 
firmed with distance estimates to the structures. There are many 
possibilities as to the nature of this overdensity. First, it could be 
part of the Sagittarius Stream, which lies in the same area of the 
Sky. However, the bulk of the Sagittarius debris lies at a lower 
declin ation (see Figure [l}. Another possibility is the Monoceros 
ring jYannv et ail 120031) . visible in the SDSS data at the same 
Galactic longitude but positive latitude. However, the kinematic 
feature in the SEGUE fields seems too localized and is l imited to 
latitudes -30° < b < -20°. Finally, dMaiewski et alj2004h reported 
the detection of an extended structure - perhaps a segment of tidal 
debris - a few degrees away from these fields as part of their survey 
of the Andromeda galaxy. Irrespective of which possibility is cor- 
rect, our hypothesis is that Segue 2 is embedded in a tidal stream. 



4 PHYSICAL PROPERTIES AND STELLAR 
POPULATION 

Assuming Segue 2 has a Gaussian velocity distribution, we mea- 
sure its velocity dispersion using a maximum-likelihood method. 
Specifically, we evaluate the marginal likelihood obtained after in- 
tegrating the usual G aussian likelihood o ver all mean velocities (see 
Equations 1 and 2 of lKlevna et aljj2~004l) ). We find the velocity dis- 
persion that maximizes this likelihood and determine the bound- 
aries that enclose 68% and 95% of the area under the likelihood 
curve. For the five bright RCGs (marked in red in Fig. [6), we mea- 
sure a velocity dispersion cr Vo = 3.4^^2 2) ' fln s '• Tms resui t ls 
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not strongly sensitive to our membership criteria - if we include the 
three fainter candidates (orange points in the Figures) passing our 
initial velocity and magnesium cuts we obtain <r Vo = 3.6*j^^2 
km s . 

Our analysis indicates that we resolve the central velocity dis- 
persion of Segue 2 , ruling out zero with more than 99% confi- 
dence. Adopting the idealised assumptions of spherical symmetry, 
dynamical equilibriu m and "mass follows light," implicit in th e for- 
mula M = 850r h cr 2 Vo jlllingworfrj[l976l : ISimon & Gehj2007h . the 
velocity dispersion of the five bright RCG members implies a dy- 
namical mass of M = 5. S^j^'^j^x 10 5 M Q , and amass-to-light ratio 

ofM/L v = 650^™ 0) [M/Lv]o. 

We can also estimate the kinematic properties of the stream 
from 15 prospective stream members satisfying -60 < V < -20 
km s _1 and 0.9 < SMg < 1.45 shown as a box in Figure [6] The 
mean velocity of the stream is -45.1 ± 0.1 (±0.2) kms~' with 1 a (2 
cr) errors. This is obtained by marginalizing over the dispersion. It 
is offset by ~ 5 kms -1 from the systemic velocity of Segue 2 . The 
velocity dispersion of the stream is 7.1^| i^^) kms -1 . 

For the brightest members of Segue 2 (the BHB and RCG 
stars), the signal-to-noise of the spectra is good enough to esti- 
mate metallicity. This was done by directly comparing the average 
continuum normalised spectrum of the three BHB spectra satisfy- 
ing 18.0 < g < 19.0 and the average of the five best RCG spec- 
tra satisfyin g 18.0 < g < 19.2, with a grid of model atmosphere 
spectra (see Wal ker et alj {2009) for further details). Using th e re- 
lationship log 10 (r eff ) = 3.877 - 0.26(g - r) i llvezic et alj|2006l) . the 
average color {{g - r) = 0.53) of the five RCG members implies 
(Trff) = 5500K for these stars, while for the BHB stars the average 
color «g - r) = -0.22) implies (T c ff > = 8600K. Anchoring the 
for the spectral fit tightens constraints on gravity and metallicity, 
giving log g = 2.5 + 0.5 and [Fe/H] = -2.0 ± 0.25 for the averaged 
BHB stars and log£ = 2.5 ± 0.5 and [Fe/H] = -2.0 + 0.25 for the 
averaged RCG stars. For the latter, each had sufficiently strong Mg 
and Fe lines and continuum signal-to-noise to model individually. 
In all cases the best fit model spectrum satisified 2.0 < log g < 3.0; 
5000K < T cff < 5500K and -2.5 <[Fe/H]< -2.0 corroborat- 
ing their categorisation as RCG stars. These colors, magnitudes, 
surface gravities and effective temperatures are all fully consistent 
with the BH B and RCG stars being at a common distance (see, e.g., 
lGravll2005l p. 57). 



5 CONCLUSIONS 

A search in the high Galactic latitude (\b\ > 20°) area covered by 
SEGUE imaging has revealed another new satellite. Segue 2 has 
a half-light radius of ~ 30 pc and an absolute magnitude of M v = 
-2.5. The photometry and spectroscopy suggest that the metallicity 
[Fe/H] ~ -2. Segue 2 is similar in structure, size, luminosity and 
velocity dispersion to three other recent discoveries, namely Segue 
1, Boo II and Coma. 

The latter three are likely to be embedded in the Sagittar- 
ius Stream. For example, Coma is superposed on the edge of the 
Stream and is at a distance that suggests association with the old 
leading (C) arm of the Sagittarius. iNiederste-Ostholt et al. I J2009h 
shows that Segue 1 stars are indistinguishable from the Sagittar- 
ius Stream, both photometrically and kinematically. It is natural to 
conclude that Seg ue 1 was once a satellite of the Sagittarius galaxy. 
iKoch et all J2009l) have shown that Boo II lies close to the young 
leading (A) and old trailing (B) arms of the Sagittarius, and has 
similar kinematics. 



There is indirect evidence that our new discovery Segue 2 is 
also immersed in a stream. To begin with, it lies on the edge of the 
Sagittarius stream, as seen by the SEGUE survey. Kinematically, 
there is a cold stream-like component in both our follow-up spec- 
troscopy and in the SEGUE spectroscopy of nearby fields. 

We speculate that all four objects are possibly satellites of 
satellites, remnants from the disruption of larger galaxies in the 
Milky Way halo. Under the simple assumption of an isothermal 
halo, the current circular velocity of Segue 2 is v c ; rc ~ 5 kms -1 . The 
original v cilc before the tiny dwarf accreted onto the Milky Way halo 
as part of a parent subgroup might have been close to 10- 15 kms -1 
(a reducton by a factor of 2-3 is expected as a result of tidal shocks 
- see Mayer et al. 2007). Initial halo masses corresponding to such 
circular velocities are below 10 7 M o . With such low masses, these 
satellites of satellites are likely to have formed before reionization, 
at z > 10, since afterwards gas would have been photoevaporated 
owing to heating by the cosmic ionizing background (Barkana & 
Loeb 1999). Such an early formation epoch would naturally imply 
a high central density and a resulting resilience to tidal disruption of 
their inner core. The inner, surviving core of the object is what we 
may be witnessing in these four objects. These, and not the more lu- 
minous, ordinary dwarf spheroidal satellites of the Milky Way and 
M3 1 , are good candidates for being fossils of the reionization epoch 
(Diemand et al. 2007). Accounting properly for their origin in the 
context of the formation of the Local Group is crucial in order to 
interpret correctly the discrepancy between the observed luminos- 
ity function of satellites and the predicted substructure mass func- 
tion. Satellites of satellites are only resolved in the most recent dark 
matter cosmological simulations of the Milky Way halo (Diemand 
et al. 2008), but not yet in more realistic fully hydrodynamical sim- 
ulations that account for the various mechanisms to which baryons 
are subject. Being able to identify satellites of satellites in cosmo- 
logical hydrodynamical simulations will thus be a crucial step for 
interpreting new observations such as those presented here. 
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